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Abstract 

By a perturbation approach, wc construct geometric solitons with various vortex struc- 
tures(vortex pairs, vortex rings) for some geometric flows(Wave maps, Shrodinger flows) from 
Minkowski spaces to C M'^. 

1 Introduction 

The well-known Landau-Lifshitz (LL) equation is a fascinating nonlinear evolution equation, 
which describes the dynamics of several important physical systems such as ferromagnets, moving 
space curves, etc. and has intimate connections with many of the well known integrable soliton 
equations, including nonlinear Schrodinger equations and modified KdV equations. It admits 
rich dynamical structures including solitons, dromions, vortices, etc. 

In many situations ferromagnetic materials may be viewed as continuous media. For instance, 
all the spins are locked together by the strong exchange interaction for a single-domain ferromag- 
netic nanoparticle at sufficiently low temperatures. In order to solve the so-called momentum 
problem in ferromagnets [22], G. Volovik proposed a simple model of delocalized electrons in [58]. 
The model introduces the normal velocity of the fermionic liquid as an additional hydrodynami- 
cal variable, describing the background fermionic vacuum, which is inspired by superfluid motion 
at T = in He^-A ([25, 26]). The purpose of this approach is the restoration of the correct linear 
momentum density by using these " hydrodynamical" variables. Moreover, for planar magnetic 
systems the existence of topologically nontrivial vortex configurations requires the fluid to be 
rotational with non- vanishing vorticity function [44] . 

To be precise, the density and the normal velocity of the fermionic liquid, introduced in [58], 
are denoted by p and v respectively. Now, since the vorticity of v is proportional to the density p 
of the topological current of the magnetic configuration, one may formulate a simple model of an 
anisotropic ferromagnetic fiuid, which arises as a modification of the phenomenological Landau- 
Lifshitz equation, with the aim of trying to solve the momentum problem. A magnetic fluid, or 
a spin-liquid, is a continuum of material particles, each of which is endowed with a magnetic 
moment. In the continuum approximation the spatial variations of the physical quantities occur 
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over a large number of particles, so that the medium is characterized by the macroscopic field 
m{r,t), i.e., the local magnetization with the constraint |mp = 1, the velocity field v(r, t) in the 
Lagrange picture, and the density of particles p{r,t). 

In a local moving frame with velocity v(r, t) = (vi,--- ,Vn), the precession motion of m = 
(mi, 7712,7713) is described by the Landau-Lifshitz equation in 77-dimensions 

9(777 + v^djUi = —777 X VE(,{m) (1-1) 

and Vi is a sort of vector gauge potential, whose tensor field, the vorticity in the hydro dynamical 
language, is given by 

diVj — djVi = 2m • {dim x djm). (1-2) 
Here Ee{S) denotes the effective energy, which is of the general form 

Ee{m) = ^ / (|V777|^ + V'/?i777j^)dx+ / m-Hdx, 

where the /3i's are the anisotropy constants, H is the external magnetic field, and the inner 
product in is defined by the constant diagonal metric g = diag(ai, • • • , a^) where aj are the 
exchange coefficients. By the definition of metric g, we have 



\{xi,- ■ ■ ,Xn)\l = '^aixj, djm = g'^ Vidjm, |Vr77|^ = ^ aj|9, 

i=l i=l i=l 



r77p. 



In the sequel, we will use the summation convection. For compressible fluid we need, in addition, 
to take account of the continuity equation 

dtp = div{pv), V = {vi,- ■ ■ ,Vn), 

and to assume that the exchange and anisotropy coefficients in Eg, depend on p. However, in the 
present paper we restrict ourselves to the incompressible case, that is /? is a constant (see [44]). 

In the case of a ferromagnetic boundary condition 

n 

m(r,t) ^ (0,0,-1) for ^ 



xf — )• 00, 



1=1 



the ri-dimensional domain space can be compactified as a topological 77-sphere S" , and 777 defines 
mappings into corresponding to homotopy classes belonging to 7r„(S^). In the case of 2 
and 3 dimensions, the mapping class [m] can be characterized by the following topological index: 

(1) for the case 77 = 2, it is the integer- valued topological degree of the mapping from onto 
(topological charge) 

(5 = — / m ■ {dim X d2m) dx; 
4vr 

(2) for the case 77 = 3, it is just the Hopf invariant of 777 from S'^ into (linking number) 



(87r)2jR3 
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where the Levi-Civita symbol for n = 3 is given by 

_ ji- -k){k- i) 

Moreover, the "magnetic vorticity" (topological current) is defined by 



Ji = -^(-iikm ■ {dj-rn X dk-m) 



which is a conserved quantity in the case of 3 dimensions and satisfies d^Ji = 0. 

In the case of one dimension the above model is just the usual Landau-Lifshitz equation[37]. 
The Ishimori model [28] can be recovered only for particular choices of the metric g = diag(l, — 1). 
It is well known that the Ishimori model is the first example of integrable classical spin model 
in 2 + 1 dimensions. 

According to the signature of the metric g we have the different cases (+,••• ,+) and 
(+,••• ,+,—,••• ,— ), possibly after a rescaling of the spatial variables. One usually denotes 
(M",5) by M-^'^ with + TV = n, and M^'^ is just the Minkowski space. In mathematics, we 
need to consider the following Landau-Lifshitz equation with hydrodynamical term (see [33]) 



dtm = my. \ {di + ■ ■ ■ + df^ - dj^^^ dj^+N)^ - Ei=i l^i^^i^i - 

-b (Ef=i Vjdjm - Ylf=K+i Vjdjm^ , 

diVj — djVi = 2m ■ {dim x djm), (1.4) 
with the incompressible condition 

div(v) = 0, 

where 6 is a real constant, ei = (1,0,0), 6*2 = (0, 1,0) and 63 = (0,0, 1). This equation unifies 
almost all Landau-Lifshitz models. 

In the present work we will be concerned only with b = and the isotropic higher dimension 
case, in the absence of an external field. Our task turn into considering the following problems 
from Minkowski space into a standard unit sphere §^ in 

dtm = m X {nm + \Dm\'^m) in M x M^'^ with m(t, r, s) G C ffi^, (1.5) 

and the equation for its stationary solutions 

Dm + \Dm\'^m = in R^'^ with m(r, s) G C M^. (1.6) 

Here, the differential operators are defined by 

am \ ^ am 

□m = Arm - Asm = - 



1 9 I 1 2 I 1 
\Dm\ = Vt-w — V,m 



dm 



dr 



N 

E 

i=i 



I dm 
Idsj 
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Note that for the convenience of future use of notation, we here denote r and s the variables in 
Minkowski space. 

It is easy to see that (1.5) and (1.6) are some types of geometric flows. More precisely (1.5) 
can be regarded as the Schrodinger flows for maps from a Minkowski space into in the sense 
of [9] (also see [53]). While (1.6) can be viewed as a generalization of the cj-model which is a 
two dimensional scalar field equation with as the target manifold, i.e., S : M^'^ — )• S^. Usually, 
one calls (1.6) for 5 : M""^'^ — )• as the wave map equation. 

From the viewpoint of analysis, many mathematicians made contributions to the Landau- 
lifshitz equation without hydro dynamical term. For the Schrodinger flows from the Euclidean 
space R" into S^, despite some serious efforts (see e.g. [1, 2], [16, 17], [20, 19], [29, 30], [21, 52]) and 
the references therein), some basic mathematical issues such as global well-posedness and global 
in time asymptotic for the equation (1.5) remain unknown. One also studied the Schrodinger 
map flow from or a compact Riemannian manifold into a Kahler manifold ([4], [9]- [11], [50] 
and the references therein) . If one is interested in one-dimensional wave (plane- wave) solutions of 
(1.5), that is, m : M X R — ^ (or S^), a lot were known as (1.5) becomes basically an integrable 
system (see [16]). The problem in 2-D or higher dimensions are much more subtle. Even though 
it is possible to obtain weak solutions of (1.5) (see [9]-[10], [13, 14], [16], [45], [56]), one does not 
know if such weak solutions are classical (smooth) or unique. 

One also studied the Cauchy problems of the Ishimori model and some existence results were 
established recently. We refer to [33] for more details. Especially, C. Kenig, G. Ponce and L. 
Vega in [32] have ever studied the following Schrodinger equation which is analogous to the 
Schrodinger flow from Lorentzian manifold: 



where u = u{t,x) is a complex valued function from R x R'^, is a non-degenerate second-order 
operator 



for some K G {I,-"" i^}; S'Hd P : C^"''"^ — t- C is a polynomial satisfying certain constraints. 
They proved the local well-posedness of the above initial value problem in appropriate Sobolev 
spaces. No doubt, it is a hard task to settle the existence and uniqueness problems of such a 
class of Cauchy problems (1.5). 

It is well-known that vortex dynamics is a natural paradigm for the field of chaotic motion 
and modern dynamical system theory. The vortex dynamics provides some physically profound 
examples of Hamiltonian systems of infinite dimensions, attracting much interest in connection 
with chaotic phenomena in dynamical systems. From the viewpoint of fluid mechanism, it is 
of important significance to study the behavior of the vortex solutions to such Ishimori type 
equation. Indeed, some exact vortex solutions for the isotropic case (/3 = 0) are considered in 
the case n > 2. For instance, in [18, 44] the two dimensional case without external magnetic 
field was studied and a wide class of solutions can be generated by using time-dependent gauge 
transformations. For the Landau-Lifshitz equation (1.1) from a standard Euclidean space R" 
with 6 = 0, /3 = and H = 0, hy our knowledge there is no existence result on vortex solutions 




— = i2u + P{u,Vu,u,Vu) 
u{Q,x) = uq{x), 



(1.7) 




j>K 
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in analytic aspects. However, one has established some existence results on vortex solutions to 
the equation with 6 = 0, /3 = (0, 0, 1) and H = Ohy employing the Lyapunov- Schmidt reduction 
method ([41, 59]). 

A natural problem is whether or not one can employ an applicable analysis method to con- 
struct some solutions to (1.5), which are of vortex structure. This is the main aim of the present 
paper. More precisely, we intend to construct some geometric soliton (or solitary wave) solutions 
with various vortex structures to the above problem (1.5) and (1.6). Thus, let's to recall some 
well-known facts and some known results on the existence of vortex solutions with finite energy 
to the Ishimori model on Ri'^ X M and the Landau-Lifshitz equation (1.3) from the Euclidean 
space with b = and the anisotropic higher dimension case, in the absence of an external 
field. 

1.1 Anti-holomorphic Ishimori model 

In the case n = 2, b = 1 and f3 = (/3i, /32, /Js) = 0, the above Landau-Lifshitz equation (1.3)-(1.4) 
without external magnetic field for m : M}'^ x R — S"^ reads 

dtm + vidim — V2d2m = m x {df — 0|)??t,, (1.8) 

diV2 — d2Vi = 2m ■ {dim x d2m), (1.9) 
with the incompressible condition 

+ ^2^2 = 0. (1.10) 

The following conservation law holds true 

Oodo - diZi + 82^2 = 0, 

where 

ao = |aim|2 + |92m|2, 
•Ui = vi2o — 2dim ■ (m x (df — 9|)m) + 2m ■ {dim x — did2m x d2m), 

and 

Z2 = 2d2m ■ {m x {df — c?|)m) + V22o — • {dfm x did2m — dim x d2m). 
The energy functional 

E{m) = / {\dim\^ + \d2m\^) dx 

and the topological charge are conserved quantities, which are related by the Bogomolny type 
inequality 

E{m) > 87r|Q|. 
It follows from the following evident inequality 

\dim lb Cijm X 9jm|^ dx > 0, 

and is saturated by time dependent spin configurations satisfying the self-duality equations 

dim ± eijm x djm = 0. 
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If we consider the standard unit sphere as the Riemann sphere for a complex plane, we have 
the following stereographic projections 



2C 1-ICI 
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mi +11712 = - r- 77T, ms 



where C,{x,y,t) is a complex value function. In stereographic projection form, the self-dual 
equation corresponding the first sign can be rewritten 

Cz{x,y,t) = dzC{x,y,t) = 

while for the second sign we have the anti-holomorphicity 

Cz{x,y,t) = dzC{x,y,t) = 0. 

Here z = x + iy and dz = ^{dx — idy), dg = ^{dx — idy). These conditions written in terms of the 
real Re^ and imaginary Im^ parts of function representing the Cauchy- Riemann equations, 
describe the incompressible and irrotational fluid flow with the velocity potential ReC and the 
stream function Im^. 

The equation (1.10) can be solved in terms of the stream function of the flow, vi = —d2f, 
V2 = diip so that we get the so called Ishimori Model 

dtm = m X (dfm — (9|m) + d2'-pdim + di(pd2Tn, 

A(p = {df + c?|)99 = 2m ■ {dim x d2m). 

In terms of complex variables, 

= vi -\- iv2 = ^iifz, V- = vi — iv2 = —2iipz 

for incompressible flow, preserving anti-holomorphicity constraint, we have dependence C, = 
C{z,t) and the model reduces to the system 

iCCf 
1 + ICP 



iCt - 2ipzCz + 2C-Z-Z - -r^^ = 0, 



We can rearrange the first equation as 

Kt + 2Cz{-^ - 21n(l + ICP) + InG-}.- = 

If we choose ip = — 21n(l + ICP), then it is easy to verify the equation (1.11) holds true automat- 
ically. Hence, we have the anti-holomorphic Schrodinger equation 

Kt + K-z-z = 0. 

By the transformation u = 4(ln(^)f, which is a complex analog of the Cole-Hopf transformation, 
it implies the following complex Burgers' equation 

iut + uuz + 2uzz = 0. 

Here, u can be interpreted as the complex velocity of effective flow with the complex potential 
f{z) = In^^. Then every zero of function (, corresponds to the vortex solution (pole of complex 
velocity) of anti-holomorphic Burgers' equation. For more details we refer to [18, 44] and reference 
therein. 
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1.2 Vortex phenomena in higher dimensional cases 

Let us consider the case of Landau-Lifshitz equation from M.^ with external magnetic fields 
(Schrodinger map flow equation) in the form 

^ = m X (Am - 771363) inMxM^, (1.12) 



or equivalently the equation 



— 777 X = A771 — 777363 + (I V?77| + m^)m. (1-13) 

Here 7/7 : R x — t- C M'^ and where 63 = (0, 0, 1) € M^. From the physical side of view, one 
expects topological solitons, which are half magnetic bubbles, exist in solutions of (1.12) (see 
[27] and [48]). Indeed, in [23]-[24], F. Hang and F. Lin have established the corresponding static 
theory for such magnetic vortices. 

By a reduction method, F. Lin and J. Wei [41] looked for a solution of the traveling wave 
777(5', Sjy — Ct) (i.e., travel in the stv— direction with the speed C = ce > 0) of the equation (1.12). 
Then m must be a solution of 

diTi 

— ce— — = m X (Am — 77^363). (1-14) 

OSN 



After a proper scaling in the space, (1.14) becomes 

-C-^ =777 X (A777 - ^63), sGM^, (1.15) 

OSN e 



or 

dm . 777363 ,9 777q, 

CmX- = A777 + ( V777r + ^)777. (1.16) 

The main results of the paper [41] are the following: Let N > 2 and e sufficiently small there is 
an axially symmetric solution 777 = 777(|s'|, sat) G C°°(M^,S^) of (1.15) such that 

/"I 2 
E,{m)= -(|V7n|2 + ^)ds < 00 (1.17) 
Jrn 2 

and that m has exactly one vortex at (|s'|,SAr) = {a^,0) of degree +1, where ^. If A'" = 2, 
the traveling velocity C ~ e, while C = {N — 2)e| log j\ for N > 3. Naturally such solution 777 
gives rise to a nontrivial (two-dimensional) traveling wave solution of (1.12) with a pair of vortex 
and antivortex which undergoes the Kelvin Motion as described above. Solutions constructed of 
this form are called traveling vortex rings for the case of the dimension N > 3. Later on, J. Wei 
and J. Yang [59] concerned the existence of traveling wave solutions possessing vortex structures 
with rotating invariance for problem (1.14). 



1.3 Main results: Special solutions with vortices for (1.5) and (1.6) 

The problems on the existence and the blow-up profiles of wave maps are well-known in analysis 
and attract much concerning in last few decades. In particular, I. Rodnianski and his coauthor 
in [49, 51] exhibit stable finite time blow up regimes for the energy critical co-rotational wave 
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map with the target in all homotopy classes and for the critical equivariant 50(4) Yang-Mills 
problem. They derive sharp asymptotics on the dynamics at blow up time and prove quantization 
of the energy focused at the singularity. For the case of degree one, a sort of different blow-up 
profiles was shown in [36]. The existence and regularity theory of wave maps from M^'^ was also 
established, for more details we refer to [55] and reference therein. However, it seems that one 
has not known whether or not (1.6) admits any solution with vortex structures by the knowledge 
of authors. 

Since we intend to seek for some special solutions of soliton type with vortex structure to (1.5) 
and (1.6), let's recall that solitons are some special exact solutions to classical non-linear field 
equations in physics. They are localized and have a finite energy. In this sense, they behave like 
ordinary particles. Solitons and multi-solitons are stable because they carry a topological charge 
N, which is an integer and equals the net number of particles. Therefore, such solitons are usually 
called as topological solitons from the viewpoint of topology. For instance, the Ishimori Model 
for the su{2) and su{l, 1) algebras can be related to the (2 -|- l)-dimensional Davey-Stewartson 
equation with pseudo-euclidean space metric [38, 40]. However, in this representation nontrivial 
mappings from space-time to inner space, classified by a topological charge, could arise. 

Because the topological charge is a conserved quantity, a single soliton cannot decay. It 
is worthy to point out that the conservation of is not due to a Noether theorem, but to 
the topological structure of the soliton. On the other hand, some geometric Hamiltonian or 
dispersive flows are of subtle symmetry when the starting manifolds and target manifolds are of 
invariance with respect to the groups of transformations. Usually these geometric flows admit 
some symmetric solutions with respect to these Lie groups of transformations, which obey some 
conservation laws due to Noether theorems and are of the properties of solitons. 

We recall that in [12] Ding and Yin proposed to study the periodic solutions of the Schrodinger 
flow in the case where the target manifold is a Kahler-Einstein manifold with positive scalar 
curvature, such a class of special solutions can be viewed as a sort of geometric solitary wave 
solution ([53]). Recently, from the viewpoint of differential geometry C. Song, X. Sun and Y. 
Wang [54] propose a notion "geometric solitons" , which is concerning a special solutions for some 
geometric flows and can be regarded as a geometric generalization of the classical solitary wave 
solutions. More precisely. 

Definition 1. [54] Suppose (M, g) is a Riemannian manifold and {N, h, J) is a Kahler manifold. 
A solution u to the Schrodinger flow 



is called a geometric soliton, if there exist two one-parameter groups of isometrics ipt '■ M ^ M , 
(j)t : N ^ N and a map v : M ^ N independent of t such that u has the form: 



For more explanations and examples for geometric solitons, the reader can refer to [54]. In 
the present paper, we will use similar ideas to transfer problem (1.5) and problem (1.6) to elliptic 
problems and then try to construct solutions with vortices. For more details, we refer the reader 
to section 2. 

Let V be the inverse of the stereographic projection 




(1.18) 



u{t,x) = (t)t o V o ipt. 



mi -|- im2 



: §2 ^ C 



1 + "7-3 
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and then set 

7^(T,•) =Ve'\ 

Here are main results. 

Theorem 1.1. [Vortex Pairs] Assume that N = 2. 

(1). For problem (1.5), there is a solution in the form 

m{t,T,si,S2) = 7^o'i/;(^,r,Sl,S2) G ^"^(IR x M x M^S^)^ 

where ip{t, r, si, S2) = U{si, S2 — cr — cot) : M x — )• C possesses two vortices of degree ±1 
traveling along the curves 

{ {xi,X2) G : xi = ±d, X2 G M}, 

under Kelvin motion with relation 

1 2c 2uj 



□ 



d 

(2). For problem (1.6), there is a solution in the form 

m{T,si,S2) = 7^o^/>(T,Sl,S2) G C°°i^ X M^S2)^ 

where iJj{t, si, S2) = U{si,S2 — cr) : M x — )• C possesses two vortices of degree ±1 
traveling along the curves 

{ {xi,X2) E : xi = ±d, X2 G M}, 
under Kelvin motion with relation 

1 2c 

d ^ 7T^' 

Theorem 1.2. [Vortex rings] Assume that N = 3. 

(1). There is a solution to (1.5) in the form 

m,{t,T,si, 82,83) = notp{t,T, 81,82,83) G C°°{R X M X W^,S^), 

where 'tp(t, r, 81,82, S3) = U{8i, 82, 83 — ct — tot) : R x M x R'^ — > C has a vortex helix directed 
along the curve in the form 

a ((icosa,dsina,0) G M^, 

and u is also invariant under the skew motion expressed by cylinder coordinate 

{r,e,S3)^ {r,e + a,83), Va G M. 

The traveling velocity and the geometric parameters of the vortex ring obeys the relation 

1 , , 4c 4a; 

: log a • 



d ^ VT 



(2). There is a solution to (1.6) in the form 

m(r,si,S2) = 7^o^(r,Sl,S2,S3) G C°°(M x M^, S^)^ 

where iP{t, si, S2, S3) = U{si, S2, S3 — ct) : R x M"^ — )• C has a vortex helix directed along the 
curve in the form 

a G M ((icosa,(isma,0) G R^, 
and u is also invariant under the skew motion expressed by cylinder coordinate 

E: {r,e,S3)^ {r,9 + a,S3), Va G M. 

The traveling velocity and the geometric parameters of the vortex ring obeys the relation 

□ 

Remark 1.3. Some words are in order to explain the results in the above theorems. 

(1) . For more information of the solutions, the reader can refer to Section 2. 

(2) . The reader can refer to Remarks 4-3 and 5.5 for the relation between traveling velocity and 

locations of vortices. In the proof these theorems, we need the traveling velocity is small 
enough 

The construction of other vortex structures such as vortex hehces is under working. The 
organization of the paper is as follows: In subsection 2.1, we introduce some transformations 
and then write problem (1.6) and problem (1.5) into elliptic cases for new unknown functions. 
To describe the vortex phenomena (vortex pairs and vortex rings), we then further formulate 
these problems by the symmetries in subsection 2.2. An outline of the strategies of proof will be 
provided in subsection 2.3. Some preliminaries are prepared in Section 3. Section 4 (Section 5) 
is devoted to the construction of solutions with vortex pairs (vortex rings), which will complete 
the proof of Theorem 1.1 (Theorem 1.2). 

2 Various vortex structures: outline of the proofs 

For the construction of various vortex structures and prove all theorems in Section 1, the main 
objective of this section is to formulate the problems in suitable forms, and then introduce some 
notation and describe the symmetries of the problems. For the convenience of readers, we also 
provide an outline of the proofs in subsection 2.3. 
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2.1 Formulations of the problems 
2.1.1 Wave maps 

Setting the stereographic projection 

mi + 11712 
1 + ms 

problem (1.6) becomes 



°^ - T^TT^ ( VrV' • - V,^ • V,^ ) = 0. 



i + IV'P 

Here and throughout the paper, we use ijj to denote the conjugate of ijj. 
We look for the traveling solitary wave solutions in the form 

V'(r, si, ■■■ , sn) = U{si, ••• , SAr)e*^, 

Sl = Si, ■■■ , SN-1 = SN-l, SN = SN - CT. 

This gives that U will satisfy 

. 1 - dU 2 d'^U 2[7 2 5C/ dU 



2/7 1 — ir/|2 

Suitable rescaling 

U{si, ■■■ , Sn) = u{si, ■■■ , Sn), 



Sl = Sl, ■■■ , SN~1 = SN-l, SN = SNVl-<^- 

will derive that 

1 — liip 2 c du . 2m ^ ^ 1 — 

1 : — TTT , — + AU — ; — rTrVli • Vti + ; — nrU 

l + \u\^ y/T^ dsN l + |n|2 l + |lt|2 

where we have denoted 

^ 9s? ' \dsi ' ' 9sAr 



In the sequel, we will set 



1 - InP 
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2.1.2 Generalized Schrodinger maps 

Setting the stereographic projection 

mi + im2 

1 + 7713 

problem (1.5) becomes 



2ip 
1 + 



idti^ + □?/'- , , r,|2 (VrV' • - V,7/; • V,^A) = 0. (2.5) 



Here and throughout the paper, we use to denote the conjugate of V'. 
We look for the traveling solitary wave solutions in the form 

^{t,T,si, ■■■ , sn) = U{si, ■■■ , SN)e''^, 

Sl = Si, ■■■ , SN-l = SN-l, SN = SN - CT - Ujt. 

This leads to the problem 

dU l-|f/P dU od^U 2U ndU dU 

— h I TTTTTT 2 C -— C ——75 h TTTTTT C 



dsN 1 + |C/|2 dsN ds% l + \U\^ dsNdsN 

2U 1-|C/P 



Suitable rescaling 



will give that 



U{si, ■■■ , Sn) = u{si, ■■■ , Sn), 

Sl = Sl, ■■■ , SN-l = SN-l, SN = snV'^ - (?■ 



bj du 1 — Itip 2 c du 
— i — , -rr: — + i ■ 



Vl - c2 dsN 1 + liiP Vl - 9sAr 

2n 1 — ItiP 

+ An r-rrrVti • Vu H = 0. 

1+up 1 + ktP 



(2.6) 



(2.7) 



2.2 Vortex phenomena 

In this subsection, for future constructions of solutions with various vortex structures, we shall 
further transform problem (2.3) and problem (2.7) to suitable forms. Note that we will assume 
that the traveling velocity parameters c and oj are small positive numbers in the present paper. 
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2.2.1 Vortex pairs in two dimensional spaces: N = 2 

We will construct solutions with pair of vortices to problem (2.3) and problem (2.7) in the case 
of iV = 2. We will set 



e = 2c/Vl-c2. (2.8) 

Problem (2.3) can be written as 

^£ ■■ , 2 + Atx + ' - — • Vw = 0. (2.9) 
We are also concerning problem (2.7) in the form 

-ZKe-— + le—^^— + + ^ ' ' n - • Vn = 0, (2.10) 

OT2 1 + |U|^ aX2 1 + |up 1 + |Mp 

where we have denoted 



Ke = a;/Vl-c2. (2.11) 

We will construct solutions with vortex pairs(a vortex and its antipair) to problem (2.9) and 
problem (2.10) in Section 4. 

Notation and symmetry properties of the operators: For later use, we will denote 

Z = Xl+ iX2 

and then set the notation 

2u 1 — \u\^ 

So[u] := An - i— ■ Vn + r-p^u, 

1 + |n|^ 1 + 

^ r 1 1 - du _ , , du (2.12) 

1 + lul"' 0x2 0x2 

Si[u] := So + eQ2, Hu] ■■= §0 + eQ2 + KffTa, 
It is easy to see the operators Si and S2 are invariant under the following two transformations 



u{z) — )• u{z), u{z) — )• u{—z). 
Thus we impose the following symmetry on the solutions u to problem (2.9) and problem (2.10) 



n := { u{z) = u{z), u{z) = u{-z) }. (2.13) 
This symmetry will play an important role in our analysis. 
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2.2.2 Vortex Rings 

Here we only consider tiie case = 3 and tlien construct solutions with vortex rings for problem 
(2.3) and problem (2.7). We set 



e| loge| 



2c 



(2.14) 



By using the cylinder coordinates 

s\ = rcos^, §2 = rsin^f, S3 = S3, 

equation (2.3) is transferred to 

, 1 — ItiP du d'^u Idu 1 d'^u d'^u 1 
z e I log e| -— — ^ — + --^ + -— + 



\u\ 



1 + 9s3 ~'~ dr'^ ~^ r dr ^ j.2 Q\^2 ~^ Qg2 ~^ ^ + 



u 



2u 



1 + \u\ 



\dr J 



\d^J \ds3J 



(2.15) 



For problem (2.15), we want to find a solution u which has a vortex ring directed along the curve 
in the form 

a G M H> ((icosa,(isina,0) G M^, 

with two parameters 



d = d/e. 

Moreover, u is also invariant under the skew motion 

S : (r, ^, S3) ^ {r,^ +0,33), V a E M. 
Hence, u satisfies the problem 

I, ,1 — |nP(9u d'^u Idu d'^u 1 — luP 

loge L , I 12 1^ + ITT + ~ir + + 1 I I 12 ' 
1 + |n|^ as3 ar^ r or as^ 1 + |u|^ 



(2.16) 



(2.17) 



2u 



(2.18) 



1 + 

For further convenience of notation, we also introduce the rescaling 

{r,h) = {xi,X2), z = xi + ix2. 

Thus (2.18) becomes 



(2.19) 



ie\ log e| 



1 — InP du d'^u 1 du d'^u 1 



1 + |n|2 dx2 ~^ dx\ ~'~ xi dxi ~'~ ~'~ 1 + 



+ 



+ 



+ 



\u\ 



\u\ 



2u 



1 + |n|2 



/ du \ ^ / du \ ^ 
\dxi) \dx2 



(2.20) 
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We will also consider problem (2.7) in the form 

II I 5ti . 1 — ItiP du . 1 — litP 2u ^ ^ 

-iK£\loge\ — h te\loge\ rio^ h An + — — p-prn = — — — njVn • Vn, (2.21) 

as3 1 + mr 0S3 1 + mr 1 + mr 



where 



Ke|loge| =a;/Vl - c2. (2.22) 



As we have done in the above, for the existence of solutions with vortex rings, problem (2.21) 
can be transformed to 

. , , du 1 — du d'^u 1 du 

-iKe\loge\ h tde\log£\-——K^ h tt^ H ^— 

dx2 1 + 0x2 ox\ xi dxi 

(2.23) 



d u 1 — \u\ 2u 



+ ^ + 



dxi 1 + ki r 1 + ki 



/ du / du \ 



\dxi ) ^ \ dx2 



The symmetric properties of the problems: Some words are in order to describe the 

methods in dealing with the existence of solutions of the above problems. Problem (2.20), and 
problem (2.23) become two-dimensional cases. The key point is then to construct a solution 
with vortices of degree +1 and their antipairs of degree —1. Additional to the computations for 
standard vortices in two dimensional case, there an extra derivative term 

\ du , 

(2.24) 



x\ dx\ 

In the following we shall mainly focus on how to deal with these terms by the methods in [41]. 

By using the symmetries, in the sequel, we shall consider problems (2.20), (2.23) with im- 
posing the boundary conditions 

\u(z)\ — > 1 as \z\ — +00, 

du (2.25) 
— (0,X2) = 0, Vx2 G M. 
dx\ 

Moreover, it is easy to see that problem (2.20) and problem (2.23) are also invariant under the 
following two transformations 



u{z)^u{z), u{z)^u{^-z). (2.26) 

Thus we impose the following symmetry on the solution u 

n := I u{z) = ^JI(i), u{z) = u{-z) |. (2.27) 

This symmetry will play an important role in our analysis. As a conclusion, if we write the 
solutions u to problem (2.20) and problem (2.23) in the form 

U(X1,X2) = Ui{xi,X2) + iU2{xi,X2), 
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then ui and U2 enjoy the followmg conditions: 

til(xi,X2) = Ui{-Xi,X2), Ul{xi,X2) = Ul(xi,-X2), 

U2{xi,X2) = U2{-Xi,X2), U2{xi, X2) = -U2ixi, -X2), (2 28) 

Note that we shall construct the approximate solution to satisfy the boundary conditions in 
(2.28). We will use these conditions to construct solutions to problem (2.20) and problem (2.23) 
in Section 5. 

Recall the notation in (2.12). We denote 



h[u]:=-^. (2.29) 

Xi OXi 



For later use, we also set 

SsM = So + e|loge|Q2 + Hi, 



(2.30) 

§4M = §0 + e|loge|Q2 + Ke|loge|T2 + Mi. 



2.3 Outline of the Proof 

To prove all theorems in Section 1, we will use the finite dimensional reduction method to find 
solutions to problems (2.9), (2.10) in Section 4, and also for problems (2.20), (2.23) in Section 5. 
The finite dimensional reduction procedure has been used in many other problems. See [6], [15], 
[39] and the references therein. M. del Pino, M. Kowalczyk and M. Musso [7] were the first to 
use this procedure to study Ginzburg-Landau equation in a bounded domain. The methods in 
the present paper basically follow those in dealing with vortex phenomena for Shrodinger map 
equations in [41]. Here are the main steps of the approach. 

Step 1: Constructions of approximate solutions 

To construct a real solution, the first step is to construct an approximate solution, denoted 
by Vd in (4.2) or (5.19), possessing a vortex located at ei = (d, 0) and its antipair located 
at 6*2 = (— d, 0). Here d is a parameter to be determined in the reduction procedure. In the 
construction of approximation to solutions to problem (2.20) and problem (2.23), there are 
singularities caused by the application of the terms in (2.24) to the phase terms of the standard 
vortices, which will be described in subsection 5.1. 

The approximate solution has the symmetry 

U2{xi,X2) = U2{xi,-X2), U2{xi, X2) = U2{-Xi, X2) ■ (2.31) 

By substituting Vd into problems (2.9), (2.10), (2.20), (2.23), we can derive the estimations of 
the errors such as §i[Vrf], S2[Vd], S4[Vd] in suitable weighted norms. The reader can refer to the 
papers [7] and [41]. 

Step 2: Finding a perturbation 
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We intend to look for solutions to problems (2.9), (2.10), (2.20), (2.23), by adding a pertur- 
bation term, say ^, to the approximation where the perturbation term is small in suitable 
norms. More precisely, for the perturbation ^ = ipi -\- i-ip2 with symmetries in (2.13) or (2.27) , 
we take the solution u in the form (cf. (4.13) or (5.25)) 

u{y) = r]{Vd + iVdij) + (1 - r])Vde'^. 

This perturbation method near the vortices was introduced in [7]. 

For given parameters d and e small, instead of considering the problems (2.9), (2.10), (2.20), 
and (2.23) we look for a tp to their projected forms in (4.31) or (5.53). By writing the projected 
problem in the form of the perturbation term ^/^(with a linear part and a nonlinear part), we can 
find the perturbation term ifj through a priori estimates and the contraction mapping theorem. 

Step 3: Adjusting the parameters 

Note that the perturbation term ip and the Lagrange multipliers c in (4.31) or (5.53) are 
functions of the parameters d. To get real solutions to problems (2.9), (2.10), (2.20), and (2.23), 
we shall choose suitable parameters d such that c is zero. It is equivalent to solve a reduced 
algebraic system for the Lagrange multiplier 

Cs{d) = 0. 

We can derive the equations in (2.32) by the standard reduction procedure. This will be done 
in subsections 4.5 or 5.5. In other words, we achieve the balance between the vortex-antivortex 
interaction and the effect of motion of vortices by adjusting the locations of the vortices. 



3 Some preliminaries 

In this section, we collect some important facts which will be used later. These include the 
asymptotic behaviors and nondegeneracy of degree one vortex. 

For the problem 

Au + ^—\^u ^^Vii-Vn = on M^, (3.1) 

1 + 1 + |np 

if we look for solutions of standard vortex of degree +1 in the form 
in polar coordinate {i,6), then p satisfies 

Another solution w~ := p{£)e~^^ will be of vortex of degree —1. These two functions will be our 
block elements for future construction of approximate solutions. 

Notation: For simplicity, from now on, we use w = p{£)e^^ to denote the degree +1 vortex. □ 
The following properties of p are proved in [23]. 
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Lemma 3.1. There hold the asymptotic behaviors: 

(1) p(0) = 0, < p{£) < 1, p' >Ofor£> 0, 

(2) p{i) = 1 - + 0(^-3/2g-£^ ^g^^ ^/^g^g Co > 0. □ 

Setting w = li^i + iw2 and z = xi + 1x2, then we need to study the following linearized 
problem of (3.1) around the standard profile w. 

Lo(0) =A(j) -I , I 12 V0 - I 12 

8('u;, (^)(wiVi(;i + 'u;2V'u;2)^ 4(Vt(;,V(^) 
+ n , I |2^2 + 1,1 12 ^ 

4{l + \Vw\'^){w,(l)) 2|Vw|2 l_|y;|2 
n I I 12^2 ^ + 111 12 '^ + 1 1 I |2 '^- 

The nondegeneracy of w is contained in the following lGmma[41]. 
Lemma 3.2. Suppose that 

LoM = 0, (3.3) 
where (p = iwtp, and -(/; = V'l + ^^2 satisfies the following decaying estimates 

l^il + l^llVV'il < C(l + |z|)-^ 

m + \z\\viJ2\ < c{i + \z\r'-<^, 

for some < g < 1. Then 



(3.4) 



, dw dw 

<P = "^1^ ^ "^2-^ — 

0x1 0x2 

for certain real constants ci, C2- □ 

Finally, we will need the following lemma on decay estimates of a linear problem in M^. 
Lemma 3.3. Let h satisfy 

Ah + f{z) = 0,h{z) = -h{z),\h\ <C (3.5) 

where f satisfies 

\m\<j^-^^^.o<e<i. (3.6) 

< (3.7) 

□ 

For the proof of the last two lemmas, the reader can refer to [41]. 



Then 



4 Vortex pairs: construction of solutions to (2.9) and (2.10) 

In this section, we will construct solutions with vortex pairs to problem (2.9) and problem (2.10), 
which will complete the proof for Theorem 1.1. 
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4.1 Approximate Solutions 

Using the degree one vortex, we introduce the approximate solutions. Throughout the paper, we 
assume that the distance d satisfies 

d=^, where [l/lOO, 100]. (4.1) 

The approximate function is then defined by 

Vd{z) ■.= w+{z-ei)w~iz-e2), (4.2) 

where ei = (d, 0), 6*2 = (— d, 0). It is easy to see that Vd^z) S IT with IT defined in (2.13). 
Furthermore, a simple computation shows that 



\Vdiz)\ 1 as \z\ +00. (4.3) 



In fact, for 1^1 >> d, we have 



{xj -cP + xl + 2ix2d) 



a/(xi - dy + xl v^(xi + dy + xl 

Here 6^ denotes the angle argument around ^ and = \z — ^\. It is easy to see that 

V^ = j^^{z-0, Ve^ = j^{z-0\ m\ = j^^, (4.4) 
where we denote z"*" = (— X2,xi). 
4.2 Error Estimates 

We plug in the approximate function (4.2) into the operators Si and §2 and obtain the errors 

El ■.= Si[Vd] = MVd] + eq2[Vd], 

(4.5) 

IE2 := HVd] = MVd] + eq2[Vd] + ^^eT2[Vd]. 
Our purpose in this subsection is to estimate these two errors. 
By our construction and the properties of p, we have 

1 + 0(e-°'^^(l^-'^il'l^-'^2|)^l_ (4^g) 



Vd = e''^^i^'^^2 



We divide M? into two regions: 

ml = {xi > 0}, ml = {xi < 0}. 

By the symmetry assumption in (2.13), we just need to consider the region 
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In the region M5_, we consider two cases. First, if — 6*2! > 2 + 



2 ' 



then we have 



Vd = w{z - ei)e"'^-2 l + 0( 



-d/2-|z-de] 



1/2) 



\Vd\ = p{\z-ei\) l + 0(e 



-d/2~\z~ei 



1/2) 



Secondly, if | z — 62 1 < | + ^^"2^ ' ^'^^ translated variable as follows 

z = ei+y. 



Then we obtain the following estimates 



VVd ■ VVd 



■ w 



l + 0{( 



Aw - 2iVw ■ VOf. 



-d/2-\y\/2\ 



e ^2 



w 



\y + ei- 6*2 



+ 0(e~^/2-H/2) 



2w 



■ e "2 X 



l + \Vd\^ " " I + 

^Q(^^-d/2-\y\/2y 

Combining the estimates above, we have for z G M^, 



-^tVw ■ V9g^ + — 



+ 0(e-"'/2-ls'l/2) 



(4.7) 



On the other hand, we can estimate the term Q2[Vrf] as follows 
dVdl-\Vd\^ 



h[Vd] = i 



dx2 l + \Vd\^ 

_d_ 

dy2 



l + 0(e 



89, 



iw- 



Similarly, 



T2[Vd] 



dw 



dy2 dy 



dVd 



£2 



1 + p^ 



dw 



e-^^?.-2 +0(e-'^/2-|y|/2). 



-d/2- 



(4.8) 



86, 



iw- 



Note that 



dy2 



dy2 dy2 
yi + 2d 



+0(e-rf/2-H/2)_ 



(4.9) 



O 



0(e), 



(4.10) 
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In summary, we have obtained for z G 



6*1 +y 



"2(p2 _ 1) 



^Vw • V6)e-2 + 2 I 1 I — ^ 



+ ie 



dw 
dy2 



ae. 



62 



dy2 



1 + 



(4.11) 



E2 = HVd] = MVd] + eQ2[Vd] + KeT2[Vd 
"2(p2 _ 1 



p2 + l 



w 



+ ^e 



dw 
dy2 



89, 



iw- 



dy2 



1 



+ l\y + ei - 62 



■ e "=2 - IKE 



e ^=2 



l + p2 



5^2 



iw- 



62 



5^2 



Similar (and almost identical) estimates also hold in the region M^. 

4.3 Setting-up of the Problem 

We look for solutions to (2.9) and (2.10) in the form 

u{y) = viVd + iVdtl;) + (1 - vWde'"^ 
where 77 is a function such that 

rj = f^{\z - ei\)+fi{\z - 6*2!). 
In the above fj is the cut-off function defined by 

r)(s) = 1 for s < 1, fj{s) = for s > 2. 



(4.12) 



(4.13) 



(4.14) 



(4.15) 



We may write 1/^ = 1^1 + ^1/^2 with "01) "02 real- valued. The symmetry imposed on u (see (2.13)) 
can be transmitted to the symmetry on tp 



(4.16) 



il;{z) = -iIj{z), iIj{z) = ipi-z). 

This symmetry will be of importance in solving the linear problems. 

By using (2.9), we can formulate the equation for the perturbation term ip. Let > 1 be 
fixed positive constant. In the inner region 



by setting 
we have 



Si = {z\z€ Bgniei) U BQR{e2) }, 

cj) = r|iVd^P + {l-v)Vd{e'^ -I), 
u = Vd + (i> 



(4.17) 
(4.18) 
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and then write the equation (2.9) in (f) as 

Li[</>] + NiM + Mi[0] = Erf, (4.19) 

where we have denoted 

The nonhnear operator is 

Ni[0] = F{Vd + 0) - F(yrf) - F'(yrf)0 + 0((1 + |0|)|V</.|2), 

(4.20) 

Mi[0] =eQ2[l"d + (/'] - eQ2[^d]. 

Note that, for problem (2.10) we have a similar formulation, but the nonlinear operator Mi is 
given by 

MiM =eQ2[Vd + 4>] - eq2[Vd] + KeT2[Vd + (P] - KeT2[Vd]. 
We now formulate the problem for the perturbation term ip in the outer region 

&2 = [z\z£ {B^Riei) U B^R{e2)y }, 
in which we have u = V^e*'^. By simple computations we obtain 



iVde'^ ^ ^d(l + |^dPe-2V'2) Vl + |yrf|2e-2V'2 

1 2|Vrfpyrf(e-^-^^ - 1) ^ 2|Vrfp(l-e-2V'2) 



i^d {l + \Vd\^){l + \Vd\^e-^^^) (l + |y^|2e-2V'2)(i + |y^|2) 

a^rf 2|yrf|2(l-e-^^^) 9V^l-^2rt^^e2^ Ei 



5x2 (l + |Kipe-2'^2)(i + |y^|2)y^ 9x2 l + |Frf|2e-2'A2 iVd' 

The above equation can be also formulated in the form 

L2M + N2M + M2[^] = E, (4.21) 

where we have denoted 

I r/1 ^ , , 2(i-|yrf|2) ^, , 4zFrfV2 . W-^^ 



l + |Frf|2 (l + |Vrf|2)2 « (l + |^dP)2' 

Vd 

^r,i . dVd 2|yrfp(l-e-^^^) l-2|Vrf|2e-2'/'2 



N2[^] = :f^V7A-VyrfO(V'2) + 0(||V^d|'-l| + |^2|)|V^-V^| + iO(|e-'^^-l + V2|), 



dx2 (1 + |Ki|2e-2V'2)(i + |y^|2) 1 + \Vd?e 



2p-2i/)2 ' 
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E 



El 



Recall that = V'l + Then setting z = ei + y, we have for z € 

/ /\i)i + o{e-\y\)\v^ 



(4.22) 



0[e~\y\\Vil) ■ VV'I + |V'2||VV • VVI + \ip2?) 



(4.23) 



/0(ei^ 




(4.24) 



Let us remark that the explicit form of all the linear and nonlinear terms will be very useful for 
later analysis. We can also write problem (2.10) as the above in a similar way, but the nonlinear 
operator M2 has the form 

dil) l-2\Vd\^e-'^^^ 



M2[ij] =2e 



2|^,|2(i_e-2'/'2) 



dVd 2|F,|2(l-e-2^^) 



+ 2ie 



2Ke 



dx2 (l + |l/d|2e-2'/'2)(i + |y,|2)y^ 



dx2 l + \Vd\^e-^^^ ' 

dip l-2|Vrf|^e-^^^ 
8x2 1 + |Vd|2e-2V'2 



2i KE 



Let us fix two small positive numbers 0<7<1, 0<£)<1. Recall that (/) = '^^^■01 
ipi + itp2 ■ Denote ij 



\z — ej\, and define 



'T(^j<2) 



(4.25) 



2 
2 



;V'i||l-(£,>2) + I|4+'V01 



7>2) 



IKj'"'V'2||L=o(,^,>2) + Pr'^^2||L-(£,>2) 



(4.26) 



tj n2\\L°°{ej>2) 

We remark that the choices of these norms are motivated by the expressions of (4.23)-(4.24). 
A direct application of (4.7), (4.11), and (4.12) yields the decay estimates for the errors 

El E2 



El 



which will be stated in the lemma. 
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Lemma 4.1. It holds that for z G (-B2_R(ei) U B2R{e2)Y 



Re(Ei), Re(E2) 



< 



Im(^i), Re(^2) 



< 



(1 + |z - ei 



+ 



[i + \z-S2\Y' 



+ 



where g G (0, 1) is a constant. Moreover, 



{l + \z-ei\y+s (i + \z-e2\y+(>-' 



(4.27) 



(4.28) 



(4.29) 



4.4 Projected Nonlinear Problem 

Let fj be defined as in (4.15) and i? > be a fixed large positive number. By defining 



dd 



R 



R 



(4.30) 



we consider the full nonlinear projected problem 

C[i)] + M[i)] + M[il:] = £ + cZd, 
Re( / 4>Zfi] =0, satisfies the symmetry (4.16), 



Ii in 61, 
II2 in ©2- 



(4.31) 



where we have denoted that 



£ = Li, M = Ni, M 
£ = L2, AT = N2, M 



Note that in the above we have used the relation (p = iVdip in ©i. It is easy to show that 
possesses the symmetries in (2.13). Here is the resolution theory. 

Proposition 4.2. There exists a constant C , depending on 7, g only such that for all e sufficiently 
small, d large, the following holds: there exists a unique solution to (4-31) and satisfies 

Ue\U < Ce'~^. (4.32) 

Furthermore, V'e,^ is continuous in d. □ 

The proof of this proposition is similar of that for Proposition 5.4. 

4.5 Reduced Problem 

From Proposition 4.2, we deduce the existence of a solution {(j),c) to (4.31). To find a real 
solution to (2.9), we shall choose suitable d such that c is zero. This can be realized by the 
standard reduction procedure in this subsection. 

On Ml, we recall that z = ei + y and 



Vd = w{y)e-'^-2 (1 + 0(e 



-^d/2~\y\/2 
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then 

Zd=^ = -^e-^^-2(l + 0(e-'^/2-H/2)) _ i^^g'^^^^^ (l + 0(e-'^/2-M/2)) , 
od oyi od ' 

where dO^^jdd is small in the sense 

In the polar coordinates, we have w = p{r)e'^^ 

— = (p'cose-i-sme)e"^, — = (p'sme + i-cos6)e"^. 
oyi r oy2 r 

Multiplying (4.31) by ^jqqp^^pp- -^d and integrating, we obtain 



Using Proposition 4.2 and the expression in (4.23), we deduce that 

Re I / ZrfNi[0] ) = o(e). (4.33) 



On the other hand, integration by parts, we have 
Let us observe that 

— So[^d] = 1 ^] = I^il^^i] = o(e), (4.34) 



dd 

and thus by Proposition 4.2 



(1 + l^dlT 

It remains to estimate the following integral 



Re( / -—^^^^U{ZM=o{e). (4.35) 



■^•^(I (TTiW "^H ^ '^(I (TTiW 
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2 Re 



+ 2eRe i 



where So[Vd] and Q2[Vrf] are defined in (4.7) and (4.8). 
The expression in (4.8) gives that 



eRe 



1 



eRe i 



eRe i 



1 — \w\'^ dw dw 
(1 + |t(;P)3 dy2 dyi 



+ o{e) 



(l + |p|2)3 



p sin 9 + i— cos ^) (p' cos 9 + i— sin 0) + o(e) 



(I-P')PP' 



+ o(e). 



r(l + p2)3 4 
On the other hand, using the estimate (4.7), we have 



Re 



Re 



Re 



1 



+1 1 + |if I |y + ei — 62] 



1^ (l + |u;|2)2 

X e"*^ (p' cos ^ + sine) 1 + o(e). 



(4.36) 



Note that in 



Hence 



Re 



r 



a + \Vd\'r 



Re 



{\w\'^-l) i 
(|u;|2 + l)3dL 



p' sin 9 + i- cos / 
r 



p cos 9 + i— sin/ 
r 



+ o{e) 



jp' - 1) P'P 

(p2 + l)3 ^ 

:(p ) + o(e) 



vr 
id 



+ 0(e) 



(4.37) 
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Combining all estimates together, we obtain the following equation 



c{d) = Co 



id ^ T 



+ o{e), 



(4.38) 



where o(e) is continuous function of d (which is a consequence of continuity of 4> in d) and 
Co 7^ 0. By simple mean- value theorem, we can find a zero of c{d). This completes the reduction 
procedure for problem (2.9), which also finishes the proof of the existence of vortex pairs of (1.6) 
in Theorem 1.1. 

For the reduction procedure of problem (2.10), we shall estimate 

/ ^ N / 



Cf-Re 



Re 



+ Re 



+ Re 



^ ZdMVd] 



1 



ZdU[4>] 



^ Zd^M+Mi 



where §2 = So + £Q2 + is defined in (2.12). The expression in (4.9) gives that 

1 



K e Re i 



K e Re i 



K e Re i 



KE 



1 dw dw 



(1 + \w\'^y dy2 dyi 



+ o{e) 



(4.39) 



(1 + |P| 



2-^2 



p sin 6 + i- cos 0) (p' cos + i- sin ^) ] + o(e 



PP 



TT 



KE - + o(e). 



r(l+/52)2 

Whence after similar computation in the above, we need consider the equation 



c{d) = Co 



vr ne e kit 



+ o{e), 



(4.40) 



where o(e) is continuous function of d (which is a consequence of continuity of (j) in d) and cq ^ 0. 
If 1 — 2k 7^ 0, by simple mean- value theorem, we can find a zero of c{d). This completes the 
reduction procedure for problem (2.10), which also finishes the proof of the existence of vortex 
pairs of (1.5) in Theorem 1.1. 

Remark 4.3. Recall the parameters e, k given in (2.8), (2.11). By the relation in (4.38), we 
have that, for problem (2.9) (i.e. (1.6)), the vortex and antivortex pair undergoes Kelvin motion 
when its speed is sufficiently small so that the vortices are widely separated with the relation 

1 2c 



d 

It is also that from (4. 40), for problem (2.10) (i.e. (1.5)), there holds 

1 2c 2uj 



d ^/l^ 
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5 Vortex helices 



In this section, we will construct solutions with vortices to problem (2.23), which will provide 
vortex rings to problem (1.5). We omit the similar arguments for the existence of solution to 
problem (2.20) (i.e. the existence of vortex rings to (1.6)). 

5.1 Approximate solutions 

5.1.1 First Approximate Solution and its Error 

Recall the vortex solutions and w~ defined in (3.2). For each fixed d := ^ with d G 
[1/100, 100], we define the first approximate solution 



vo{z) := w'^ {z - ei)w (2-62), 
where ei = {d,0) and 6*2 = {—d,0). A simple computation shows that 

|uo(-2)| — )■ 1 as \z\ — )■ +00. 



In fact, for |z| >> d, we have 



gif <?! - if 62 



(xf -d'^ + xl + 2ix2d) 



^J{xi - d)2 + xl ^J{xl + dY + xl 



^, d 
1 + 0^ + ^1. 



|zp \z\ 

Here 9^ denotes the angle argument around ^ and = \z — ^\. It is easy to see that 



where we denote z-^ = {—X2,xi). 



1 



By our construction and the properties of p, we have 



I _|_ (9(e-min{k-eiMz-e2|)^ 



We divide the region & into two parts: 

6+ = |^:xi>o}, ©_ = |z:xi<o}. 



(5.1) 



(5.2) 



(5.3) 



(5.4) 



(5.5) 



By our symmetry assumptions, we just need to consider the region S+. In the region 6+, we 
consider two cases. Firstly, if |z - 62! > f + then we have 



^^0 

l^ol 



w{z - ei)e-'^^'2 1 + 0(e"rf/2-k-'iei|/2) 
p(|z-ei|)fl + 0(e^'^/2-l"-"^il/2; 
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Secondly, if |z — 6*2! < ^ + -^^2^' ^^^^ ^^^^ paper, we often use the translated variable 

as follows 



z = ei+y. 

The term IHIi['yo] obeys the following asymptotic behavior 



(5.6) 



xi dxi 



d 



d + yi dyi 



w 



1 



yi 



where we have 



d d{d + yi) 
d . 



iw- 



dw 



dyi dyi 



-X2 



dyi (y, + 2d)^ + yl 



0(d-2) = o{e^). 



Whence, by (5.3), the term ™ ^1 ^ singularity in the form 



/(«^0) ~ - jT|2- 

d dyi d\y\^ 



(5.7) 



(5.8) 



In summary, by recalling (4.7), (4.8), and (4.9), we have obtained for z G (5+ with z = ei+y 
§4N = SoN] + e|loge| Q2N] + Ke|loge| T2[wo] + IHIi[t;o] 



2{p'^ - 1) /j^ _ 1 ^ 



+ l 

+ ie\ loge|e~*^'='2 
— iKe\ log e\e~^^^2 



p2 + 1 |y + ei - 6*2! 



50, 



dw 



dy2 dy2 



iw- 



+ 0(e-'^/2-|2;|/2) 



|!!i_i^^ + 0(e-'^/2"M/2) 
9^2 9y2 



(5.9) 



The similar (and almost identical) estimates also hold in the region z € (3_. A direct application 
of (5.9) yields the following decay estimates for the error 



SqN] + £|log£| Q2N] + K£|logg| T2M 



< 



Lemma 5.1. It holds that for z G (i?2(ei) U i?2(e2))' 

Re(£^4 

Im(^4) < 
where g £ (0,1) is a constant. 



+ 



(l + |z-ei|)3 (1 + 1^-621)3^ 



+ 



{l + \z-ei\y+e (i + |^_e2|)i+f^ 



(5.10) 

(5.11) 
(5.12) 
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Proof: The estimates for the term SoboJ^:^ follows from (5.9). We just need to estimate the 



terms 



e|loge| Q2M- 



1 

ivo 



Ke\ loge| T2['yo]- 



whose definition are given in (2.12) and (2.29). 
Let us compute for z G S+ 

0X2 I 0X2 fo 



dp{\z-ei\) 1 , . d{es,-e. 

£ ZTTT + 



dx2 p{\z - 6*1 



dx2 
xi — d 



xi + d 



(xi - dy + xl {xi + dy + xl 



The estimate (5.11) then follows. Let us notice that for z G (3+, \z — ei| < d, 



xi — d 



Xl + d 



< Ce- 



1 



{xi-dy + x^ {xi + dy + x^ (l + |z-ei 



On the other hand, for z G S+, \z — ei\ > d, we then have 



Xl — d 



Xl + d 



{xi — dy + X2 {xi + dy + x'2 



< c 



< c 



{l + \z-ei\y+^' 



{l + \z-ei\y - (l + |z-ei|)i+f' 



Thus (5.12) is proved. 



□ 



5.1.2 Further Improvement of the Approximation 

As we promised in previous subsection, we now define a new correction ip^ in the phase term. 
The phase function ipd will be decomposed into two parts: singular part and regular part. Then 
we will define an improved approximation and estimate its error by substituting to problem 
(2.23). 

To cancel the singularities in (5.8), we want to find a function $(2/1,2/2) by solving the problem 
in the translated coordinates (yi,y2) 



+ 



in 



dyf dyj d7^|yp 
In fact, we can solve this problem by separation of variables and then obtain 

1 



^(yi,y2) 



4dT 



-y2 log|y| 



(5.13) 



(5.14) 



Let X be a smooth cut-ofi' function in a way such that x{z) = 1 for z G B_d^{ei) and x = 

10 

for z G (Bd{ei)) . Let (pd{z) = ^s{z) + ^r{z)-, where the singular part is defined by 



t \ ( \ 1 k - ei 



4d-f^ 



\Z - 62 



2 ■ 



(5.15) 
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While by recalling the operators Hi in (2.29), we find the regular part ^Priz) by solving the 
problem 



A 



A + : 



Note that the function (ps is continuous but Vips is not. The singularity of (ps comes from its 
derivatives. 

By simple computations, we see that for z G B_d^{ei), 

10 



A + 



4x2 (xi - d) 



0{d-^) = 0{e^). 



+ 



X2{xf - X2 - d^) 



XI ((:ci-d)2+x2) ((xi + d)2+x2) 



For z G {B_d_{ei))^ , it is easy to see that we also get 0{e'^). In fact for z G (i?d(ei))^, (^9^ = 

10 5 

and 



A + Mi 



-4X2 



((Xi - d)2 + x2) ((xi + (i)2 + x2) ' 

Going back to the original variable (r, S3) in (2.19) and letting ip{r, S3) = (pr{z) we see that 

C 



Ar,s^,^ + Mil 



< 



(5.16) 



Thus we can choose (fr such that (3 = O ( 1 ) . The regular term ojj. is in the original 

variable (r, S3). 

We observe also that by our definition, the function 

:=%-% + + V'r, (5.17) 

satisfies 



A + 



ip> = on 



(5.18) 



From the decomposition of ipa, we see that the singular term contains X2 log 1 5;— ei | which becomes 
dominant when we calculate the speed. 



Finally, we define an improved approximation 

Vd{z) :=w+{z-ei)w-iz-e2)e''^''. 

5.2 Error Estimates 

Let 



(5.19) 



p = p{\z - ei\)p{\z - e2\), <p 
So Vd = pe'**^. Note that for xi > 0, we have 



-'ei "e2 



p{\z-e2\) = 1 + 0(6 



rl^-eil 



)• 



(5.20) 
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Since the error between 1 and p{\z — e2\) is exponentially small, we may ignore p{\z — e2\) in the 
computations below. We shall check that is a good approximate solution in the sense that it 
satisfy the conditions in (2.28) and has a small error. 



Let us start to compute the errors: 



Here we have used the fact (5.18). We continue to compute other terms: 



Wd ■ VVd 



|Vp|2 - p2|V(^p + 2ipVp ■ V(p 



We then obtain that 



-ipSo[(p] + 



p2 + 



^(|V^|2-|V0p) 



+ 2i- ^Vp- + O e-2 — 

1 + p^ V 



In a small neighborhood of ei, we write 



J_dVd 

xi dxi 

Xl X\ 

ipni(p + 0{e). 



The estimates for the terms 



Ke| loge|T2[Vd] = iK£\\og£ 
e\\oge\Q2[Vd] = i£\\oge 



dVd 
dx2 ' 
1-W dVd 
1 + W 9X2 



are same as before. 
The total error is 



= ^^f— ^ (IV^P - |V0|2) + 2ii-^Vp . + 0(e 



+ ze| log £\ 
+ i£\ log el 



2 2 



dy2 dy2 



iw- ^ 



dy2 dy2 
Here we have used the relation in (5.18). 



1-/0 

e-*^--2 + 0(e|loge|e'"'/2-l^l/2) 



(5.21) 



(5.22) 



(5.23) 
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Setting z = ei + y, we then have 



Vips = -^^logdVy2 + Oielog\y\), Vy., = 0(e). (5.24) 



Thus 



|V(^|2-|v0p = o{\Vips-ve\) + \Vip,\'^ + o{e\ve\), 

Vp-Vtf = O {e\ log e\p') + 0{ep log \y\). 
These asymptotic expression wih play an important role in the reduction part. 

5.3 Setting up of the Problem 

Now we introduce the set-up of the reduction procedure. We look for solutions to (2.20) and 
problem (2.23) with boundary condition in (2.25) in the form 

u{z) = rj{Vd + iVdi^) + {l-ri) Vd e'^ , (5.25) 

where is a function such that 

r? = r?(|z-ei|)+ry(|z-e2|), (5.26) 

and fj{s) = 1 for s < 1 and r){s) = for s > 2. 

The conditions imposed on u in (2.25) and (2.27) can be transmitted to the symmetry on ip 

^(z) = i^i-z), i^iz) = -W), ^(0'^2) = 0. (5.27) 
More precisely, for ■0 = ■01 + iip2 , there hold the conditions 

'lpl{xi,X2) = 1pl{-Xi,X2), '(pl{xi,X2) = -lpl{xi,-X2), 

1p2ixi,X2) = 1p2{-Xl,X2), 1p2ixi,X2) = 1p2ixi, -X2), 



(5.28) 



This symmetry will be important in solving the linear problems in that it excludes all but one 
kernel. We may write ip = ipi + 11^2 with 1^1,1^2 real- valued and then set 

u = Vd + (l), ^ = r]iVdi; + (1 - v)Vd{e'^ - 1). (5.29) 

In the sequel, we will derive the explicit local form of the equation for the perturbation term ^|;. 
Let -R > 1 be a fixed constant. In the inner region 

(3i = {z\z£ BQR{ei) U ^9^(62) }, 

we have 

u = Vd + (t>, (5.30) 

and the equation for <j) becomes 

U[4>] + Ni[0] = Ed. (5.31) 
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In the above, we have denoted Unear operator by 

urn = A,^ + - -ij^VKi • V0 - -^^VVd ■ VVd 

+ (i + |y,|2)2 • VF, + F (F,),^. 

The nonhnear operator and the error are 

NiM = F{Vd + </.)- F{Vd) - F'{Vd)(t> + + |0|)|V</.|2 

, ■ XII I , • XII I ^'^ 

+ «Ko loge e— — + zKo loge £7- — , 
0x0 0x2 



(5.32) 



and 

Erf = -S4[Vd]. (5.33) 
In the above, we have use the definition of F in (2.4). 
In the outer region 

62 = { z I 2 G (S4H,(ei) U BiR{e2)y }, 
we have u = V^e^^. By simple computations we obtain 

+..|log.|..^^ 



^ ■ All I 1 - iVdpe-^^^ Erf 
+ iKd\ loge e- ——5 — 5:^—7^ h . 

We can also write the problem as an equation of = V^i + ^V'2 

+ N2['0] + M2['0] = E, (5.34) 
with conditions in (5.28). In the above, we have denoted 

. 4|yrf|2^2 



= :f^VV'-VyrfO(V'2) + 0(||Frf|2-l| + |V^2|)|V^-V^| + iO{\e-^^ -1 + IIJ2\) 

, • All I 

+ z loge e- — . 

0x2 
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There error term has the form 

iVd 

Recall that "0 = "01 + *'02- Then setting z = ei + y, we have for z S 

AV'i + iri& + oie-\y\)m 
+ irSf - iM^^^ + o(e-i^i)v^2 J 



(5.35) 



N2[^] 



0(e-ls^l|V^- V^l + |^2||V^-VV| + |t/'2p 



N2[^] 



(5.36) 



(5.37) 



Let us remark that the explicit form of all the linear and nonlinear terms will be very useful for 
later analysis in resolution theory. 



Let 



E 

E4 = S^iVd], 1E4 = 



(5.38) 



Based on the form of the errors, we need to use suitable norms. Let us fix two positive numbers 

p> 13, < Q <l. 

Recall that (p = iVdil^, ip = ipi + iip2- Denote £1 = \z — ei\ and £2 = \z — e2\, and define 



= \\iVdh\ 



LP(6i) + X] [ l|-^j'^^^llli°°(S2) + I|4'^^^2||l<-(62) 



In the above, 



VK2.p(6i) + Yl 



+ 



V'2||l-{62) + ll^j'^^W2||L-{62) 



©1 = {zG6:|z — ei|<3or|z — e2|<3}, 
&2 = { z e 6 : \z - ei\ > 2 and \z - 6*2! > 2}. 



(5.39) 



We remark that we use the norm L^^^ (or W^^^) in the inner part due to the fact that the error 
term contains terms like elog |z — ei| which is not L°°-bounded. 

Using the norms defined above, we can have the following error estimates. 
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Lemma 5.2. It holds that for z G &2 



|lm(E4 



< 



+ \z- eil 



+ 



+ 



[l + \z-e2\r' 



{l + \z-ei\Y+s (1 + 1^-621)1+^?' 



and also 



where g G (0, 1) is a constant. As a consequence, there holds 



5.4 Projected Linear and Nonlinear Problems 

Let 



(5.40) 
(5.41) 
(5.42) 

(5.43) 
□ 



and the co-kernel 



+ TT^ + 



dxi xi dxi 

d\ 



+ 



2(i-|y,i2. 



l + |V^d|^ 



VVd -VVd- i 



dd 



R 



R 



(5.44) 



where fj is defined at (5.26) before. Then satisfies the symmetry (2.27). 



As the first step of finite dimensional reduction, we need to consider the following linear 
problem 



-^o(^) = h in &, ip satisfies the symmetry (5.27), 



Re 



(5.45) 



6 



dx) =0 ior (j) = iVdip in Bi(ei) U Bi{e2). 



We have the following a priori estimates. 

Lemma 5.3. There exists a constant C , depending on g only such that for all e sufficiently 
small, d ~ ^, and any solution of (5.45), it holds 



* < c\\h\u,. 



(5.46) 



Proof: The proof is similar as in Lemma 5.1 in [41]. Suppose that there exists a sequence of 
e = e„ — )• 0, functions /i„ which satisfy (5.45) with 



^ II* — 1) ll^nll** — ^(^)' 
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We will derive a contradiction by careful analysis of the estimates. 

We derive inner estimates first. We have the symmetries and boundary conditions for ipi 
and V'2 in (5.28). Whence we may just need to consider the region 

= { xi > 0, < X2 < A7r/(e7) }. 

Then we have 

Re(^J^JnZd^ =2Re^^ ^nZd^ =0. 

Let z E z = ei + y and cpniu) = 4>n{z)- Then as n — )• +oo, 

Vd = w+{y)e-''^2[l + 0{e~'"^)) = -w+{y)+o{l). 
Since HV'"!!* = !> we may take a limit so that 4>n. — 4>q in Mfg^, where 4)^ satisfies 

Lo[(/>o] = 

where Lq is defined by (3.3). Observe that 4>q satisfies the decay estimate (3.4) because of our 
assumption on tp"^. By Theorem 3.2, we have 

dw dw 
90 = ci-^ h C2^. 

oyi oy2 



Observe that <^o inherits the symmetries of (j) and hence 0o = (t>Q{z). (The other symmetry is 

dw 
dy2 

dyi ' 



not preserved under the transformation z = ei + y.) But certainly does not enjoy the above 
symmetry. Hence (po = ci^^. On the other hand, taking a limit of the orthogonality condition 



we obtain 



Re I 1^ cPnZdj = 0, 



Rel / 00^1 = 0. 



This implies that ci = and hence we have 

(t>n - 

which implies that for any fixed R > 0, 

2 



inEL 



^(^\\<Pi\\Lp{ej<R) + \\4'2\\Lp(ej<R) + \\^4>i\\Lp(e,<R) + l|V(/>2||LP(^^.<i?)) = o(l). (5.47) 

We use the L^-estimates in the inner part {jz — ei| < R}. By choosing p > N large we obtain 
the embedding W^^^ into C^^'^ for any a G (0, 1). 

Next we shall derive outer estimates. Let f/ be a cut-off function such that fj{s) = 1 for s < 1 
and fj{s) = for s > 2. We consider the new function 

V' = ^'Xiz), where x(z) = 1 - v(}—^^) - ^(^~4~^)- 
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Using the explicit forms of Lq in (5.35), the first equation becomes 
^ + + ^ = 0(e-l^l)|V^| + O(VxVV'i) + O(ViAx) + h.x- (5.48) 

On the region S\(i?4(ei) U -64(62)), we have the conditions in (5.28). Moreover we have 

iVx-VV-l = o(l)(|z-ei|2 + |z-e2|2)-^. 
For the outer part estimates, we use the following new barrier function 

Biz): = Bi{z) + B2{z), 

where we have denoted 

Bi{z) = 4x'^ + 4x'^, B2{z) = Ci{l + \z\')-'^\ 
In the above /3 -\- 1/ = —g, < g < v < 1. 
Now, we do the computations for Bi 



^^^J_dB2 ^ _cc,(l + \z\')-'-'/\ 
dxl dx'2 xi dxi 

where Ci depends only on (3 and v. On the other hand, there holds 



1 dBr ^ 



l^^~^{xi-d)+l^f^{xi + d) 



xi dxi x\ 
Thus for i\ < Cgd or £2 < Cgd, where Cg is small, we have 

1 dBl ^ ^ I n-2-e . n-2-g' 



< C2C,{i^'-' + 



xi dxi 

where C2 depends only on f3 and v. For ii > Cgd, there holds 

-^<c(i+\z\'r'~'/\ 

xi 0x1 

By choosing Ci large, we finally have 

d'^B ^ ^ J_dB_ ^ _^.^_2_g^^_2-^,N 
dxf dxi ^1 ~ 2 ; • 

By comparison principle on the set Il\(i?4(ei) U -B4(— ei)), we get that 

\i^i\<CBi\\h\U, + o{l)), VzG S\(54(ei)US4(-ei)). (5.49) 

Elliptic estimates then give 

2 

EII4^''I^^iI|Il-(£,>4) ^ cm** + oil)). (5.50) 
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To estimate ip2, we perform the same cut-off and now the second equation becomes 
dxl xi dxi dxl (1 + |Vd|2)2 

= 0(-^)VV'i + 0(e~l^l)VV2 + 0(VxVV2) + 0(AV2) + /i2X- 

1 + \y\ 

Note that we also have the conditions in (5.28). Since for z G T\[Bi{ei) U-B4(— ei)), there holds 



+ -4' 
by standard elliptic estimates we have 

U2\\l^(i^>a) < C{U2\\l^(i^=^){1 + + \z- eil + \z- dei\)-^-s, 

(5.52) 

|V^2| < C(||V2||loo(,^,=^))(1 + ||V'||*)||/i||**(l + N - eil + \z- dei|)-2-^. 

Combining both inner and outer estimates in (5.47), (5.49)-(5.50) and (5.52), we obtain that 
= o(l), which is a contradiction. □ 

We consider the full nonlinear projected problem 

C[^] + + = £ + cZd, 
Red ^Z,)=0, ^satisfies the symmetry (5.28), 



(5.53) 



where we have denoted that 



£ = Li, 7^ = Ni, M = Ml in ©i, 
£ = L2, = N2, M = M2 in ©2- 

Note that in the above we have used the relation (j) = iVdtp in &i. By contraction mapping 
theorem, we make a conclusion by the following the resolution theory 

Proposition 5.4. There exists a constant C , depending onp, q only such that for alle sufficiently 
small, d large, the following holds: there exists a unique solution ipi^^ii to (5.53) and ipir^d satisfies 

\mU < Ce'-^. (5.54) 

Furthermore, 1^^,(1 is continuous in d. □ 

Here, we do not give the proof to the last proposition. The reader can refer to the arguments 
in [41] for details. 
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5.5 Reduced Problem 

From Proposition 5.4, we deduce the existence of a solution {ip,c) to (5.53). To find a real 
solution to problem (2.23), we shall choose suitable d such that c is zero. This can be realized 
by the standard reduction procedure in this subsection. 

Multiplying (5.53) by p^qp^^pp-^d and integrating, we obtain 



''''' ' (i + m?^'"' 



Using Proposition 5.4 and the expression in (5.32), we deduce that 

Re(^^Z^NiM) =o{e). 
On the other hand, integration by parts, we have 



(i + |y,|2)2 ^L^v vy^, (i + |y,|2)2 



Let us observe that 

-dVr. 



d_ 

dd 

and thus by Proposition 5.4 



As the strategy in standard reduction method, we are left to estimate the following integral 
-I^ef / -—^—^Z~dEA=Re([ , I, .... Z^^S^Va]) (5.55) 



The expressions of these error terms are given in (5.23). We will estimate the above integrals in 
the sequel. 

On (3+, recall that z = ei + y. Then we have the estimates for the first term in (5.55) 



Je+ iX + P ) \d + yidyi p^ + l ' J oyi 

f 1 2p(l -p2) ^ dd ^, ^ 
Je+ (1 + P ) 1 + r dyi 
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Let us notice that 



Hence 



2d 



log dVy2 + 0{elog\y\). 



1 , dp 



0{e) 



and also 



pip^ - 1) / ^ ,^ ,2N dp 2^(1 - " 



1 



log d 



P 

l-p^ pp' 



+0{e) 



d " Je^{p'' + lf\y\ 



dy + 0{e) 



8d 



\ogd + 0(e). 



Recall the estimates in (4.36) and (4.39). the last terms can be estimated by 



Ke\ log e\ Re 



^ ZdT2[Vd] 



-Ke\ loge|-, 



e| loge| Re 

As a conclusion, there holds 
c{d) = cqtt 



+ (l + |V^dP)2 



el log el 



1 1 — 2k 

— log d + — ^ — e| log e| + 0{e) 



(5.56) 



where cq 7^ 0. Therefore, for any 2k — 1 < 0, we obtain a solution to c{d) = with the following 
asymptotic behavior: 



— logd~2(l-2K)e|loge| 



(5.57) 



Remark 5.5. Recall the parameters e, k, given in (2.14), (2.22). By the relation in (5.56), we 
have that, for problem (1.5), the vortex ring travels along S3 axis when its speed is sufficiently 
small with relation between it geometric parameters and traveling velocity 



- log d 
d 



4c 



Similarly, for problem (1.6), there holds 

1 
d 



log d • 



4c 
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